A Magnetic Model of the Tetragonal-Orthorhombic Transition in the Cuprates 
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It is shown that a quasi two dimensional (layered) Heisenberg antiferromagnet with fully frustrated 
interplane couplings {e.g. on a body-centered tetragonal lattice) generically exhibits two thermal 
phase transitions with lowering temperature - an upper transition at Tto ("order from disorder 
without order") in which the lattice point-group symmetry is spontaneously broken, and a lower 
Neel transition at Tn at which spin-rotation symmetry is broken. Although this is the same sequence 
of transitions observed in La2Cu04, in the Heisenberg model (without additional lattice degrees of 
freedom) (Tto ~ Tm)/Tm is much smaller than is observed. The model may apply to the bilayer 
cuprate La2CaCu06, in which the transitions are nearly coincident. 
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I. INTRODUCTION 

In addressing the physics of high temperature super- 
conductivity (HTC), an important, but often overlooked, 
issue is the relation between the lattice structure and the 
electronic physics. In the highly studied "214" family of 
HTC superconductors, there is a structural phase tran- 
sition from a high temperature tetragonal (HTT) to a 
low temperature orthorhombic (LTO) phase with a tran- 
sition temperature, Tto(x) which drops as a function 
of doped hole concentration, x, from Tto(O) ~ 500K in 
undoped La 2 Cu04, to Tto(0.22) = in La 2 _ a; Sr. l; Cu04 
with x > 0.22 0,0, II Hi- We were initially motivated 
to examine the problem of interlayer magnetic coupling 
because, in the Neel state, the coupling between near- 
neighbor planes is frustrated in both the T-phase of the 
214 compounds [e.g., LSCO and Sr 2 Cu02Cl 2 0), and 
the T'-phase (e.g. Pr 2 Cu0 4 or Nd 2 Cu0 4 ) @) due to the 
location of the Cu 2+ ions on a body-centered tetragonal 
lattice. The HTT- LTO transition of the T-phase removes 
the frustration, and always occurs well above any mag- 
netic ordering temperature, Tn- On the other hand, no 
lattice distortion has been identified in the T' systems, 
even though the transition temperature to the magnet- 
ically ordered Neel state is similar for compounds with 
the two structures 0, • 

The question we address is how a purely electronic sys- 
tem, with no lattice coupling, could spontaneously break 
the point group symmetry which results in the perfect 
cancellation of the magnetic coupling between the sub- 
lattices. This approach is contrary to conventional wis- 
dom, which holds that the HTT-LTO transition in LCO 
is driven by a lattice parameter mismatch between the 
Cu0 2 planes and the interstitial layer La-O. There are, 
however, reasons to question whether this explanation 
is complete: 1) Doping with Ba ions could be expected 
to relieve the stress more efficiently than doping with 
Sr, because of its larger ionic radius. It doesn't seem 
to do that: the Tto — > at roughly x=0.20 in both 
cases. 2) The in-plane magnetic susceptibility is found 



to be spectacularly anisotropic 0, |n| below Tto {Xbb 
is between two and three times larger than Xaa, where a 
and b refer to the orthorhombic a and b axes); were the 
electronic anisotropy simply a response to the symmetry 
breaking lattice distortion, then especially at tempera- 
tures above Tn but below Tto, one would expect the 
anisotropy of any electronic response to be small in pro- 
portion to the magnitude of the orthorhombic distortion, 
which is around 4% in these materials. 



We consider undoped La 2 Cu04 where, due to the 
strong repulsions between electrons, the low energy elec- 
tronic physics is well approximated by a quantum Heisen- 
berg antiferromagnet with a single spin 1/2 on each Cu 
site. We define a continuum field theory (a O(N) non- 
linear sigma model with N = 3) in terms of the cor- 
responding long- wave- length excitations, and solve it in 
the large N limit (i.e. in the self-consistent phonon ap- 
proximation) . We show that the spin model undergoes a 
similar sequence of two phase transitions upon cooling as 
does La 2 Cu04 - a high temperature T-0 transition and a 
lower temperature transition, Tn < Tto, a t which spin- 
rotation invariance is spontaneously broken resulting in 
a Neel state. (This corresponds to the schematic phase 
diagram in Fig. 2c.) The T-0 transition in this model is 
an Ising transition and is an example of "order by disor- 
der without order" , as we explain below. The magnetic 
order is collinear, consistent with what is determined for 
La 2 Cu04. However, the difference between the two tran- 
sition temperatures is very small, (Tto — Tn)/Tn <C 1, 
while what is observed in La 2 Cu04 is a substantial dif- 
ference, (Tto — Tn)/Tn ~ 1; thus, it is clear that the 
lattice degrees of freedom play a very significant (possi- 
bly dominant) role in the stabilization of the LTO phase 
in these materials. On the other hand, Tto an d T y are 
very close in the bilayer cuprate La 2 CaCu06 an d 
the model may apply in this case. 
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II. THE MICROSCOPIC MODEL 

We start by defining a model of M coupled layers, each 
of which consists of a square lattice Heisenberg antifer- 
romagnet, with the sites on one layer situated above the 
plaquette center of the layer below: 



M 



H = EEE^K«- S 



(i) 



2 n,H u,R+5q 
11=1 R 5 
M-l 

+ E E E Jn S n,R ' S n+1,R+8 1 
71=1 R 5i 

where S n ^ is the spin at site R in plane n and the 

sums over So and Si run, respectively, over pairs nearest- 
neighbor sites within a plane and in neighboring planes: 
So = ±x, ±y (in units in which the in-plane lattice con- 
stant = 1) and 5i = ±(1/2) [x + y], ±(1/2) [x - y}. Be- 
low, we will consider the bilayer (M = 1) and the body- 
centered tetragonal (M — > oo) versions of this problem. 
We will also consider cases in which the spacing between 
the layers are not all the same so that J' n is n dependent, 
but we will always assume J J' > 0. 

In a classically Neel ordered state, S n ^ = M„ expjiQ- 

R} with Q — (7r, 7r), there is no coupling between the 
staggered order parameters, Mi and M 2 , in the two 
planes. However, high energy short wave-length quan- 
tum fluctuations induce an additional effective interac- 
tion between the spins at low energy, H — ► H + Hind- 



Hi, 



M-l 

n=i a 



n,R 



(2) 



where K,r = Eii [e lQ ' 5l S^- S n+1 ^ +5 -J is a composite 
order parameter field and [28j K n ~ (J' n ) 2 / J . For K n > 
0, the staggered moments on neighboring planes tend to 
be colinear to maximize the expectation value of $. 

Rotation by it/ 2 about an axis through any lattice site 
R on plane n is a symmetry of the Hamiltonian. Such 
a rotation transforms (f> ^ — > — $ ^, so any state with 
($) is a state of broken symmetry. For instance, the 
state pictured in Fig. 1 has ($ 1 s) < 0, but upon rotation 
by 7r/2 about a site in the lower plane, one obtains the 
state in which all the spins in the upper layer are reversed, 
and hence (& 1 changes sign. More generally, a non- 
zero ($ 1 fi) occurs only in a state in which the discrete 
C4 rotational symmetry of the lattice is broken down to 
C2, i.e. in an "Ising Nematic" state, or equivalently an 
electronically orthorhombic state. 



III. ORDER FROM DISORDER 

In two dimensions: A mechanism is called "order by 
disorder" when it is driven by fluctuation induced terms 



that are absent at the mean-field level 

ei 13 m A clas- 
sic example is the bilayer of Heisenberg antiferromagnet 
shown in Fig. l(the case M = 1). Here, in addition to 
the usual ground-state degeneracy which is implied by 
the spontaneous breaking of spin-rotation symmetry, at 
mean-field level there is an added (continuous) degener- 
acy due to the frustrated character of the intersystem 
coupling; the orientation of the staggered magnetization 
on the upper layer can be rotated by an arbitrary angle 
relative to that on the lower layer at no cost in energy. 
However, quantum fluctuations (Eq. |2J lift this degener- 
acy and produce a preferred relative orientation such that 
the staggered magnetization on the two planes is colin- 
ear. As we have seen, this state also breaks the discrete 
rotational symmetry of the lattice. 

The fluctuations that lead to order by disorder are pre- 
dominantly short-wavelength, high energy quantum fluc- 
tuations, and hence they produce the same tendency for 
the spins on neighboring planes to be colinear, even when 
there is no long-range spin order. For instance, since the 
bilayer problem is a 2D problem, at any finite tempera- 
ture, T, spin rotation symmetry cannot be spontaneously 
broken. However, the Isin g Ne matic order can (and does) 
survive to finite T 

[UHli. Thus, although originally 
conceived as being a consequence of small fluctuations 
about a magnetically ordered state, the tendency of these 
fluctuations to produce this order persists even when 
long-wavelength fluctuations have destroyed the antifer- 
romagnetic order, itself, leaving £s(T) < 00 -hence E3 
"ordering by disorder without order." 

We can estimate the transition temperature as being 
the point at which the induced interaction between cor- 
related blocks of spins is of order T (in units in which 
k B = 1) 



T, 



(bi) 



TO 



(3) 



exp[2AirJ/T], the ordering temperature, 



Tin ~ J/log[J/fC\ is significantly smaller than J only 



Since £s 

,(6iJ 

if K is astronomically small. (Here A w 1/4 for spin 
1/2 due to quantum corrections.) The resulting phase 
diagram is shown in the schematic phase diagram (a) in 
Fig. El 

In three dimensions: The magnetic Hamiltonian on 
the HTT lattice is the three dimensional extension of the 
same model, where M — > 00 and J' n and K n are inden- 
pendent of the layer index, n. However, to place the 
problem in a broader context, in Fig. 2 we have consid- 
ered a somewhat more general version of this Hamilto- 
nian in which J ^> J' 2n +\ = J' > J^n = J" > and 
K n = K(J n /J') 2 . For J" I J' = 1, this corresponds to 
the LTT structure of La 2 Cu0 4 , while for J" /J' = this 
corresponds to an array of decoupled bilayers. 

The behavior of the phase boundaries for 1 > J" / J' > 
can be readily deduced from the single bilayer results by 
scaling. The Neel temperature rises steeply from zero as 
T N - 2 At:. J I log [J/ J"}. The T-0 transition rises from its 



finite bilayer value as Tto ~ Tr 



(bi) 



TO 



J{J"/J) 2 /~<, where 
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7 = 7/4 is the susceptibility exponent of the 2D Ising 
model. 

To obtain a general solution valid for J' <~ J", we can 
no longer perturb about the isolated bilayer, and so must 
rely on other, approximate methods. 



IV. CONTINUUM O(N) MODEL 



We start with the following continuum Hamiltonian, 



J 



I' I' M 1 

Hc= dxdyJ2h(Ps(\^M 2 + |W>„| 2 ) - gi {<j> n • V>„) 2 + y&AfyVn - 92(<Pn ■ 4>n-l? + y^A^n-l] (4) 

n=l 



r 



where ip n {r) and <j>n(f) are the local staggered magne- 
tization order parameters in even and odd layers and 
\<f> n {x, y)\ 2 = \ip n (x,y)\ 2 = 1. Here the spin stiffness 

p s cx J, gi cx ^j-, (72 oc ^-j-, Vi ^' and ^2 oc J". (We 
have considered multiple coupled bilayers, but by mak- 
ing g and 77 dependent on layer index, the more general 
model can be studied as well.) We also take the order pa- 
rameters to be N component vectors, generalizing from 
the physical N = 3, since in the large N limit (which we 



expect to be qualitatively correct for N = 3) the phase 
diagram of above model can be obtained as follows. The 
partition function is given by 



Z = J D\ l D\ 2 DaD(j)Dijexp{- j J dxdy^L n ) (5) 



where the Lagrangian is given by 



~N 



XKl^nl 2 - 1) + A 2 (|^„| 2 - 1) + 2[a ln (<j> n • V„) + a 2n ((p„ ■ Vn-i)] 

+Tll4>ndxdyll)n + m&ndxdylpn-l] 

I 



T< + ^ + ^[p s (|va 



9i 
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l 2 + |vvg 2 ) 

(6) 



where \ l n {f),i = 1,2 are the Lagrangian multiplies for 

ip n {r) and <j) n (r) respectively and (Ji ni i — 1,2 are the 
Hubbard-Stratonovich field. The saddle point of above 
Lagrangian is determined the following self-consistent 
equations by taking <7j„ (r) = crj , i = 1,2 and A„ (r) = A 
(subject to the stability condition 2A > \a\ + cr 2 1): 



Ci = 



(27T) 
1 = TO 2 + 



, , d 2 kdk z (cosk z Goo — (— l) i smfc z Goi)(7) 
T 



(2tt 



1 3 



d z kdk z G 



(8) 



where we have taken A^(r) = A, i = 1, 2, m is the stag- 
gered moment which is non-zero only in the magnetically 
ordered phase, the integral for k has an ultra-violet cut- 
toff A, and 



G(k x , k y , k z 



A(k) B(k)\ 
-B(k) A(k) ) 



(9) 



where A(k) = X+ + =(<t+ +r]+k x k y ) andB(fc) = 
^ s%nkz_ ((j — +f]^k x k y ) with cr ± = (Ti ±cr 2 and ?? ± = 771 ±772 



First, let us look the uniform case g\ = 52 = g, 771 = 
r)2 = i] and o\ = a 2 = o. If 77 = 0, the above self 



consistent equation can be solved exactly. There is a first 
order transition at T c given by 



9 — 2PsA -4npJTr 
4tt Ps T c c 

For T > T c , these equations have only one solution: m = 



a = 0,A = 



-4ir Ps /T 



of solutions with m 
T 

non-zero magnetization, 



< T c , the minimal solution has a 



At T = T c , there is a family 
and A + VA 2 - a 2 = For 
A = a/Ann* and 



m 



(T c - T) _ T£_ 
T c Anp s 



log[T/T c 



(10) 



The transition at T c is a peculiar: The orthorhombic 
order parameter jumps to a finite expectation value, 
$ cx a = g/AT{p Sl just below T c , making it a first order 
transition. At the same time, the magnetization grows 
continuously from zero as T drops below T c . (Unsur- 
prisingly, the coincidence of the two transitions is lifted, 
as we shall see now, by non-zero 77.) Of course, because 
A = a, there is a gapless (Goldstone) mode for all T < T c . 

When 77 7^ 0, there are two phase transitions. The first 
one is a second order Ising-typc transition at Tto and the 
second one is the Neel transition at Tjv. Tto > ?iv for 
any nonzero value of 77. The Ising transition is marked 
at Tto where a (Tto) = and ^\t=t to 7^ and the 
Neel transition happens at Tn where A — a. From these 
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conditions, equations which determine Tto and Tn can 
be derived. While it can be proved that the solutions of 
the equations exist and satisfy Tto > Tn when 77 ^ 0, it 
would be tedious to present the proof. 

However, when r\ is small, we can solve the critical 
temperatures up to the second order of 77. Expanding 
the self consistent equations up to the second order of 77 

exp(- 

and defining a = 4^-, /? = and j(T) = 

the two critical temperatures are determined by 



TO) 



3/3 2 -/(Tn)- 1 (3 2 , . 

lo a 6 



From these equations, it is transparent that Tto > Tn- 
Up to the second order in /3, the difference between Tto 
and Tjy is given by 



5T = Tto — Tn 



48(4^p s - T N ) 



(12) 



While this result shows that the fluctuations will in- 
duce the Ising order transition before the temperature 
reaches the Neel transition, 

As 77 is expected to be small, it follows that for the 
uniform case, so is (Tto~Tn)/Tn- For nonuniform cou- 
pling parameters 171 ^ (72 and rji ^= 772, there are also two 
phase transitions and the difference between two transi- 
tion temperatures can be large. In the large N limit, 
the Ising transition critical temperature is mainly de- 
termined by <7i and rji while the Neel transition criti- 
cal temperature is mainly determined by 52 and 772. To 
show the difference between Tto and Tjy in this dimer- 
ized model, we take J = 0.12ev, r)\/p s = 0.1, p s = JS 2 , 
and 9i = 0.26J5A^|. We choose the cutoff A = 100 
and S — 1/2. By varying 772, we find that the transition 
temperature Tto — 470fc, which is almost independent 
of 772- In figGl we plot Tr °~ Tjv as a function of ln^. 
It fits to a very good linear dependence. So far, there is 
no experimental evidence supporting layer dimerization 
along c-axis in cuprates. The quantitative results here 
can not be directly linked to LCO. However, it will be 
interesting to see whether the results can be applied in 
other systems. 

If just for a test of the accuracy of the large N method 
used here, and as a comparison to the result in two di- 
mension, it is interesting to apply the same method to 
the bilayer problem. The continuum Hamiltonian for this 
case is H c in Eq. 01 with M = 1. The corresponding 
saddle-point equations are obtained from Eqs. [7| by re- 
placing the integral over k z with a sum over k z = 0, 
7T. These equations are solved with m — at all non- 
zero temperatures. There is, however, a non-zero criti- 
cal temperature for the T-0 transition, Tto, such that 
a = for T > T TO , and A = (p s A/T) exp[-4n A^pjT] 
where A„ = -Jl — (r)T/2p s ) 2 . Taking this expression for 
A(T), the critical temperature is the solution of the im- 
plicit equation Ai:A ri p s X = g, from which it follows that 



Tto ~ 47r j o s /log[(7/p s A]. As T is decreased below Tto, 
<j rises continuously from 0. 



V. DISCUSSION 

In La2- x Sr x CuOi : the HTT to LTO transition oc- 
curs at roughly twice Tn, and so is inconsistent with 
our findings. Probably, this signifies a significant role 
of the lattice degrees of freedom in driving the transi- 
tion. However, this does not necessarily imply that the 
electronic considerations found here are entirely unim- 
portant - in particular , th e order one anisotropy in the 
magnetic susceptibility [TjJ in the temperature range be- 
tween Tn and Tto argues in favor of a significant elec- 
tronic contribution to the transition. A more compelling 
case for an electronically driven transition to an or- 
thorhombic phase can be made for the two-layer com- 
pound La2CaCu2C>6, where the magnetic and structural 
transitions are observed and they are close together in 
temperature 0,0]. 

It is a very general feature of the fluctuation induced 
coupling between frustrated planes that K > 0, and 
hence that the induced order is collinear. However, a 
finite concentration of vacancies can induce the opposite 
sign of interplane coupling which can therefore lead to 
non-collinear ordering of magnetic moments in different 
ravers|l5j. When a vacancy is created in one plane, the 
LTO order near the vacancy is destroyed. The effect of 
vacancies on Tto is illustrated in Fig. 0] At mean-field 
level, the net field applied on any spin in one plane due 
to the spins in the plane below is zero; however, where 
there is a vacancy in the lower plane, all the neighboring 
spins in the next plane feel a net field which is J' times 
minus the spin that was removed to make the vacancy. 
As is well known, an applied field in an antiferromagnet 
favors the ordered moments in the plan perpendicular to 
the field, which permits a net canting of the spins in the 
direction of the applied field - this is the driving force for 
the usual spin-flop transition. Thus, there is a negative 
contribution to the average coupling, K vac ~ —n vac J' 2 /J 
which is proportional to the vacancy concentration, n vac . 
We note the AF systems i? 2 Cu0 4 (R= Nd,Pr,Gd,Eu) 
IH El E3 forming in the T structure exhibit the 
combination of no detectable orthorhombic distortion 
and non-collinear ground states. Possibly, the difference 
of the T' compounds from La2Cu04 is the moment as- 
sociated with the rare earth atoms ( for Eu atoms, due 
to the relatively small energy gap between the first ex- 
cited state and the ground state, they carry significant 
magnetic moments at 300K although the ground state is 
nonmagnetic) ; including them leads to a magnetic Hamil- 
tonian different from Eq. (1). 
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FIG. 1: A sketch for the bilayer Heisenberg model in square lattice. 




FIG. 2: a) The phase diagram for a bilayer system: Tm = Tto ^ 0. b)The phase diagram for uniform case with r\ — 0. The 
transition is first order and T to = Tn c) The general phase transition. There are two second order transitions: Tto > Tjv- 
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FIG. 3: The difference between two critical temperatures Tto T Tjv as a function of ln(gi/g2) 
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FIG. 4: A sketch of spin configurations in presence of one vacancy. When a vacancy is created in one plane, spins in the 
nearest neighbor plane favor non-collinear alignment. If alignment is collinear, which is illustrated by dashed arrows, the third 
and fourth spins are not happy energetically. 



